How good are the Garvey-Kelson predictions of nuclear masses? 
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The Garvey-Kelson relations are used in an iterative process to predict nuclear masses in the 
neighborhood of nuclei with measured masses. Average errors in the predicted masses for the first 
"q^" ■ three iteration shells are smaller than those obtained with the best nuclear mass models. Their 

quality is comparable with the Audi- Wapstra extrapolations, offering a simple and reproducible 
procedure for short range mass predictions. A systematic study of the way the error grows as a 
, function of the iteration and the distance to the known masses region, shows that a correlation 

exists between the error and the residual neutron-proton interaction, produced mainly by the 
implicit assumption that V np varies smoothly along the nuclear landscape. 
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I. INTRODUCTION. 

Mass is the most basic property of atomic nuclei. Knowledge of nuclear masses is essential not only for nuclear 
physics but for other branches of physics. In astrophysics they are required as a basic input for the calculation 
of diverse processes in the evolution of stars Several methods and models have been developed to predict 
the unknown nuclear masses Q, all of them trying to use the essential physical information determining the 
experimentally known masses. Unfortunately the predictions of these approaches tend to diverge from each other, so 
there is a permanent search for better theoretical models that reduce the difference with the experimental masses 
and produce reliable predictions for unstable nuclei. 

^| ■ Besides the global formulas, there are a number of local mass formulas. These local methods are usually effective 
, when calculating the mass of a nucleus, or a set of nuclei, which are fairly close to nuclei of known mass, taking 
f^S ' advantage of the relative smoothness of the masses as a function of proton and neutron numbers to deduce systematic 
■ * | trends. Among these methods, there are a set of algebraic relations for nuclear neighbors known as the Garvey-Kelson 
relations ||. 

o\ : 

The surface in the neutron number N, proton number Z space of the known nuclear masses, as well as the surface 
of its "derivatives" : the two-neutron and two proton separation energies, and the Q Q and Qp values, vary smoothly 
. rH , as a function of N and Z. This property is fundamental and permits extrapolation to the unknown region using only 
the systematic trends of this surface [4| . This method provides the best short range mass extrapolations Q , which 
have been published for sets of three or four nuclides in the neighborhood of those with measured masses in the 
Atomic Mass Evaluations 5, 6]. These predictions are performed nucleus by nucleus, combining a (rather elaborate) 



graphical analysis with relevant physical information 
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In the present article it is shown that the well known Garvey-Kelson relations can be used in an iterative way to 
predict masses as accurately as in AME, and significantly better than any nuclear mass model, for the first three 
iteration shells, which cover a large number of nuclei. After a short review of the Garvey-Kelson relations, compar- 
isons with several mass models are presented, as well as a study of the systematic behavior of the error in this relations. 

The Garvey-Kelson relations @ are algebraic expressions that relate the masses of six neighboring nuclides. These 
relations can be derived in an (extreme) single-particle shell model by constructing a suitable combination of neigh- 
boring masses in such a way that the residual neutron-neutron, proton-proton and neutron-proton interactions cancel 
out. The main assumption is that the positions of the single particle levels and the residual interactions between 
nucleons vary smoothly with atomic number. Two algebraic equations can be obtained in this way 

M(N + 2,Z-2)- M(N, Z) + M(N, Z-l)- M(N + l,Z-2) + M(N + 1,Z)- M(N + 2, Z - 1) = (1) 
M (N + 2,Z)- M(N, Z - 2) + M(N + 1, Z - 2) - M (N + 2,Z-1) + M(N, Z-l)- M (N + 1, Z) = (2) 



where M(N, Z) is the nuclear mass. None of the nuclei in a Garvey-Kelson relation should have N — Z odd, and for 
relation [7] to be fulfilled it is necessary that N > Z, to secure a proper cancelation of the isospin dependence of the 
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FIG. 1: Number of Garvey-Kelson estimates used to predict each nucleus. 



residual interactions [|J . 

For each relation it is possible to obtain the mass of a given nucleus in terms of the remaining five masses involved, 
which can be done in six different ways. Using the two relations it is possible to obtain up to 12 estimates for the 
mass of a nucleus @, Q . To estimate a mass using these relations it is necessary to know at least five of the six masses 
involved [1, Q. Although only a rather small set of nuclei with unknown masses satisfy this condition, it is possible 
to make further predictions through an iterative procedure, by employing the masses computed in previous steps. In 
this article we analize in a precise fashion how the average error behaves as a function of the number of iterations. We 
show that it is possible to achieve very reliable mass predictions for the first three iteration steps, which add up to a 
large number of nuclei. 



II. SHORT RANGE PROBES. 

In order to compare the predictive power of the Garvey-Kelson process with the Audi-Wapstra extrapolations for 
the AME95[5( compilation, we start with the experimentally known masses in AME95 as the initial set, predict the 
new masses reported in the AME03 compilation Q, and compare them with their measured values. 

The Garvey-Kelson relations have been used before in an iterative way to predict masses of nuclei not too far 
from known ones. A study of their predictive power for this test have been published[2]. However, to the best of our 
knowledge, an iterative process considering the 12 possible estimates for each mass has not been analyzed. Figure 
[5] shows the number of estimates used to predict each nucleus for the test described before. While most of these 
nuclei are predicted with only one estimate, it is important to note that it is just one of twelve possibilities. We also 
observe that many nuclei are predicted with more than one estimate. Given the simplicity of the method, we be- 
lieve that these two factors make it worthwhile to make a new analysis of predictability for the Garvey-Kelson relations. 

At each iteration step we estimate the mass for all the nuclei in the prediction set for which it is possible to use at 
least one Garvey-Kelson relation. If it is possible to use more than one relation, we take the average value of these 
estimates. For all the masses computed in this fashion we calculate four relevant quantities at each iteration, the root 
mean square deviation (a 1 and u % mod ) for the nuclei predicted on iteration i. 
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where the index j goes over the iVj nuclei predicted on iteration i, ei_j is the experimental error and Wj is defined 
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where the mod versions of a and Ao~ are the model standard deviation ■ The predicted masses are then added to 
the "known" set for the next iteration. 

Tableland Figure^ allow a comparison, at each iteration, between the a 1 and Aa % mass deviation for Garvey-Kelson 
and Audi-Wapstra. 

In the second column of Table [J the number of nuclei whose masses are predicted up to each iteration is given. 
Between parenthesis is the number of nuclei predicted at each iteration. From the 389 new masses, 93 of them are 
predicted in the first iteration, and 198 in the first three iterations. We see that for the first three iterations the 
Garvey-Kelson relations provide mass predictions which are comparable to those computed in Ref. Q. The quality 
of these predictions is significantly better than those provided by the best global models [2j. 

It is interesting to compare this predictions with the previous ones[2j. In the cited reference they mention "in the 
limit of seven iterations, for example, 242 out of the 382 new masses could be calculated with deviation Aa mo d = 0.232 
MeV. A total of 21 iterations were required to reach 340 of these masses, but with the resulting large rms deviation 
of Aa mo d — 0.717 MeV". Following the iterative process used in this paper, 242 nuclei are predicted after only 
5 iterations and are predicted with an Aa l mod about 0.35 MeV. To reach 340 predicted nuclei only 11 iterations 
are needed and the deviation with the experimental values is Ao l mod — 0.628. The inclusion of the twelve possible 
estimates allows the iterative process to reach further on each iteration and because of this the error grows more slowly. 

The Garvey-Kelson iterative process is a simple, systematic and reliable option to estimate nuclear masses which 
doesn't require an elaborate graphical analysis, in the vicinity of the experimentally known region it is competitive 
with Audi-Wapstra extrapolations. After the fourth iteration the iterative process tends to increase rapidly and its 
predictions are no longer reliable. The color code in Figure [3a] shows the nuclei predicted at each iteration step for 
this test. It would be useful to develop new strategies to allow the Garvey-Kelson relations to reach further. 

Fourteen different tests have been developed in order to gauge the predictive power of different mass models 10]. 
One of them, the 95-03 probe, is particularly well suited to test short range predictions and will be used to compare 
the ability of the Garvey-Kelson relations to predict masses generated by some of representative global models. 
Starting with the AME03 data set Q, we divide it into two subsets, one assumed to be known, while the other is 
predicted. We analyze the cases with N > 8, Z > 8 and V > 28, Z > 28. 

The Garvey-Kelson predictions do not vary if we consider N>8,Z>8orN> 28, Z > 28 because of its local 
character. However, some of the mass models improve their predictions when we do not include the light nuclei. In 
other words, the a % curve as a function of distance oscillates less if the light nuclei are not present in the fit. 

For the 95-03 test with N > 8, Z > 8 there are 1760 nuclei in the initial (fitting) subset, while the prediction 
consists of 389 nuclei. The case with N > 28, Z > 28 has in the initial subset 1454 nuclei and 371 in the prediction set. 
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FIG. 2: Comparison between the a 1 (a) and Aa % (b) deviation in MeV for the predictions of the nuclei on each 
iteration for the Garvey-Kelson iterative process (black) and the Audi-Wapstra extrapolations (blue) as a function 

of the accumulated number of nuclei with predicted masses. 
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FIG. 3: Set of nuclei whose masses are predicted on each iteration for the AME95-03 probe. 



The performance of the Garvey-Kelson relations (GK) to predict masses will be measured against the following 
three mass models 10]: the Bethe-Weizsacker Liquid Drop Model (LDM), modified to include volume and surface 
contributions in the symmetry energy and a Wigner correction [ll|, [T§]; the same LDM including schematic shell 
corrections, linear and quadratic in the number of valence nucleons (LDMM) [l2|, [H| and the Duflo-Zuker (DZ) 
model [IH with 33 parameters. The LDMM and the DZ models contain information about the single-particle 
(shell) behavior of nuclei through microscopic corrections, while the LDM is purely macroscopic. In each model the 
parameters arc fitted to the nuclei in the initial set for each probe. The model then produce a prediction, which 
can be compared to the corresponding GK prediction. The comparison with the Duflo-Zuker model is especially 
interesting since it produces the lowest errors against measured masses in g lobal calculations. 

Figure [4] shows the a 1 and Aa % deviations in MeV for the predicted nuclei on each iteration for the three models 
in the case N > 8, Z > 8. It is evident that the Garvey-Kelson procedure is at least as good as the models, including 
microscopic corrections, LDMM and DZ, up to iteration 7 and is better than the LDM up to iteration 13. We can 
also observe that the predictions of the microscopic models are quite stable through all the iterations, especially 
for the Duflo-Zucker model. Figure [5] shows the a 1 and Aa l comparison for the case with N > 28, Z > 28. Again 
the Garvey-Kelson procedure improves the results up to 7 and 13 iterations, respectively, when comparing with the 
precision of the LDMM and DZ models and with the LDM, respectively. 
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FIG. 4: a % (a) and Act 1 (b) deviation (in MeV) for the mass predictions on each iteration, for the AME95-03 probe 
with N > 8, Z > 8, of the Garvey-Kelson iterative procedure (black), LDM(blue dotted), LDMM (green dotted), 

DZ(red dotted). 
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FIG. 5: The same as Fig. Hfor N > 28, Z > 28 . 



III. LONG RANGE PROBES. 

It has been demonstrated previously and confirmed systematically in this work, that the Garvey-Kelson relations 
are an excellent tool to predict nuclear masses in the vicinity of the experimentally known region. However, the rms 
deviation grows rapidly because of the use of previously predicted masses on each new iteration. It is important 
to develop strategies to make reliable predictions further away from known masses. We now present a systematic 
analysis of the GK errors and their patterns of increase. 

To analyze the error produced on each iteration the AME95-03 probe is not optimal, We require a test covering a 
larger extent of nuclei, as it is very important to control the position of the predicted nuclei on each iteration. This 
ensures that the iterations truly correspond to distance to the known region. Each of the predictions of previously 
used models (LDM, LDMM and DZ) is a good set of data consisting of approximately 9000 binding energies that 
satisfy the Garvey-Kelson relations to a good approximation. Using these sets we construct the long range probe, 
again dividing the data into two subsets. All nuclei with A > 160 and (N — Z) < 44 are used as known masses and 
are shown in figure [6] in blue. Nuclei with A > 160 and (N — Z) > 44 are the ones that will be predicted by the 
iterative procedure and are shown in red. 

Figure [7a| shows the nuclei predicted on each iteration for this test. Note that in each iteration only nuclei 
with N-Z=constant are predicted. In this way, the distance to the known region and the number of iterations 
are the same and it is possible to analyze how the error grows as a function of this two variables. We cal- 
culate the accumulated root mean square deviation Act 1 between the different models and the Garvey-Kelson 
prediction, for each iteration. Figure [7b] shows Act 1 for these three models. We observe that for each of the three 
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FIG. 6: Nuclei used as input (Blue) and those who were predicted (Red) using the Garvey-Kelson relations iterative 

procedure in the long range probe. 



sets of data used, the rms grows systematically, the Duflo-Zucker model being the one for which the error grows faster. 

Figures [Sal Hi and rrUal show the distribution of the error, M model (N, Z) - M GK (N, Z), for LDM, LDMM and DZ 
sets of data respectively. It is remarkable that these errors are systematic and well localized. Taking into account 
that the three models satisfy the Garvey-Kelson relations remarkably well, we may assume that the systematic error 
is largely due to the iterative procedure. 

One of the main assumptions of the Gravey-Kelson relations is that the residual interaction between nucleons vary 
smoothly with N + Z and N — Z. In order for the different masses to cancel, the Garvey-Kelson relations must be 
closely related to the interaction energy of the last i neutrons and the last j protons in a nucleus 

e ln - JP (N + j, Z + i) = B(N + i,Z + j)- B(N, Z + j) - B(N + i,Z) + B(N, Z) (8) 

Especially for ei„_i p , ei„_2 P , £in-ip and 62 n -2p, it is even possible to derive the Garvey-Kelson relations from 
assumptions about this quantities [15j. In fact, it can be shown that one of the relations assumes that ein-ip is 
independent of N + Z and the other assumes that ei„_i p is independent of N — Z. It has been shown[l|| that this 
assumption is not followed by the experimental masses and it produces systematic errors in the predictions obtained 
with the Garvey-Kelson relations. If ei„_i p does not vary smoothly, then on each iteration of the Garvey-Kelson 
process the residual interactions do not cancel completely and a systematic error accumulates. Similar behavior has 
been observed previously in the average interaction of the last protons with the last neutrons SVn V = \tin-2p and it 
has been shown that this quantity not only is not smooth but e xhibits very distinctive patterns [17]. 

Figures [8bl I9bl and 1 1 Obi show 5V np calculated from data of each model. The LDM shows a very smooth variation of 
5V n p with N + Z and N — Z, this is the reason why the error of the iterative process grows very slowly. Comparing 
the error and SV np we can observe that both present the same trend and display some correlation. The variation of 
SV np in the LDMM is not so smooth, it has discontinuities at the shell closures and near the middle of them. This 
discontinuities are produced because the model has terms that take into account the number of valence nucleons and 
in this regions the valence particles change nature from "particles" to "holes" and vice versa. Figure [Sal shows that 
this discontinuities, where SV np varies abruptly, are boundaries for the regions which display systematic errors. 

The comparison between the errors and the 6V np in the DZ case shows more clearly the influence on the errors 
produced by the sudden variations of 5V np . The DZ model introduces a competition between two behaviors, for 
each shell, in order to model the dependence of mass on deformation. The boundary between these two regimes 
creates a discontinuity that is enhanced in the SV np . Figure [lObl shows clearly one of this lines as well as the same 
discontinuities due to valence nucleons in the shell closures and in the midshells. Both kinds of discontinuities create 
boundaries for the regions with systematic errors. 

It is clear that there exists a correlation between the regions with systematic error and the residual neutron-proton 
interaction. At least two effects produce this correlation, the dependence of SV np with N + Z and N — Z and the 
discontinuities or abrupt variations in 6V np . The discontinuities form boundaries and the error correlates with them. 
These discontinuities have been observed in the behavior of the models and it is clear that the errors of the iterative 
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process are correlated with their position. Figure [TT] and the analysis in [TtJ shows that the experimental data also 
presents discontinuities in 5V np . It is an important task to include these effects in order to improve the GK procedure. 




FIG. 7: (a) Set of nuclei whose masses are predicted in each iteration for the long range test, (b) Aa l deviation (in 
MeV) for the mass prediction in the long range test using the masses predicted by different models, LDM(blue 

dotted), LDMM (green dotted) and DZ(black). 




FIG. 8: (a) Errors for the Garvey-Kelson iterative process using the predictions of LDM. (b) Average interaction of 

the last protons with the last neutrons 5V np . for LDM. 



IV. CONCLUSIONS. 

In conclusion, a systematic procedure to estimate nuclear masses close to the experimentally known region has 
been presented. This procedure uses the Garvey-Kelson relations through an iterative process. The predictive power 
of this procedure is measured as a function of the number of iterations in the process. 

The predictions of the Garvey-Kelson iterative process have been compared with the predictions of several models: 
the Audi-Wapstra extrapolations, the Bethe-Weizsacker Liquid Drop Model, the Bethe-Weizsacker Liquid Drop 
Model with schematic shell corrections in the number of valence nucleons and the Duflo-Zuker model, as a function 
of the distance to the known region. Special attention was given to the Audi-Wapstra extrapolation technique since, 
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FIG. 9: (a) Errors for the Garvey-Kelson iterative process using the predictions of LDMM. (b) Average interaction 

of the last protons with the last neutrons SV np . for LDMM. 
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FIG. 10: (a) Errors for the Garvey-Kelson iterative process using the predictions of DZ. (b) Average interaction of 

the last protons with the last neutrons SV np . for DZ 



like the GK relations, it is based on the smoothness and the systematic trends of the mass surface in the N-Z plane, 
and because this method has the best estimates for nuclear masses in the short range. 

In all cases the Garvey-Kelson mass predictions are better than those obtained from any other model for the first 
three or four iterations. This procedure seems particularly well suited to predict unknown masses of very unstable 
nuclei, starting from a set which can include the more than 400 nuclei whose masses were recently measured at the 
nuclear laboratories at GSI, ISOLDE, GANIL, MSU, Jyvaskyla, Stockholm and FSU OH. 

A systematic study of the way the error grows as a function of the iteration and the distance to the known masses 
region shows that a correlation exists between the error and the residual neutron-proton interaction, produced mainly 
by the implicit assumption in the Garvey-Kelson relations that V np varies smoothly along the nuclear landscape. The 
discontinuities on V np produce systematic errors that grow with the number of iterations. 

While the use of Garvey-Kelson relations is one of the simplest and more accurate methods to predict masses near 
the experimentally known region, it is necessary to develop strategies to increase the range of reliable predictions. 
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FIG. 11: Average interaction of the last protons with the last neutrons SV np for AME03 experimental data. 



These strategies must take into account the fact that V np does not always vary as smoothly as the Garvey-Kelson 
relations assume and thus it is necessary to design ways to correct this effect. 

Thanks are due to C.Thibault for very useful comments. Partial financial support from Conacyt-Mexico and 
DGAPA-UNAM is acknowledged. 

TABLE I: a 1 and Aa % deviation in MeV for the predictions on each iteration of the Garvey-Kelson iterative 
procedure. The first column shows the iteration number, in the second column the number of predicted nuclei up to 
that iteration, in parenthesis the nuclei predicted only on that iteration. A graphical comparison of the data is 

shown in Fig. [5] 
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0.4012 
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